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Abstract 

A binary matrix can be scanned by moving a fixed rectangular window (sub- 
matrix) across it, rather like examining it closely under a microscope. With each 
viewing, a convenient measurement is the number of Is visible in the window, which 
might be thought of as the luminosity of the window. The rectangular scan of the 
binary matrix is then the collection of these luminosities presented in matrix form. 
We show that, at least in the technical case of a smooth mx n binary matrix, it can 
be reconstructed from its rectangular scan in polynomial time in the parameters m 
and n, where the degree of the polynomial depends on the size of the window of 
inspection. For an arbitrary binary matrix, we then extend this result by determining 
the entries in its rectangular scan that preclude the smoothness of the matrix. 

Keywords: Discrete Tomography, Reconstruction algorithm. Computational com- 
plexity. Projection, Rectangular scan. 



1 Introduction and Definitions 



The aim of discrete tomography is the retrieval of geometrical information 
about a physical structure, regarded as a finite set of points in the integer 
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square lattice Z x Z, from measurements, generically known as projections, 
of the number of atoms in the structure that he on hues with fixed scopes 
(see [5] for a survey) . A common simphfication is to represent a finite physical 
structure as a binary matrix, where an entry is 1 or according as an atom 
is present or absent in the structure at the corresponding point of the lattice. 
The challenge is then to reconstruct key features of the structure from a small 
number of scans of projections [7], eventually using some a priori information 
as convexity [1] [2], and periodicity [3]. 

Our interest here, following [6], is to probe the structure, not with lines of 
fixed scope, but with their natural two dimensional analogue, rectangles of 
fixed scope, much as we might examine a specimen under a microscope or 
magnifying glass. For each position of our rectangular probe, we count the 
number of visible atoms, or, in the simplified binary matrix version of the 
problem, the number of 1 in the prescribed rectangular window, which we 
term its luminosity. In the matrix version of the problem, these measurements 
can themselves be organized in matrix form, called the rectangular scan of the 
original matrix. Our first objective is then to furnish a strategy to reconstruct 
the original matrix from its rectangular scan. In the sequel, we will address 
this problem to as Reconstruction{A,p,q), where A is the rectangular scan, 
and p and q are the dimension of the rectangular windows. As we also note, 
our investigation is closely related to results on tiling by translation in the 
integer square lattice discussed in [6]. 

To be more precise, let M be an m x n integer matrix, and, for fixed p and 
q, with 1 < p < m,l < q < n, consider a p x q window Rp^q allowing us to 
view the intersection of any p consecutive rows and q consecutive columns of 
M. Then, the number Rp^q[M)[i, j] on view when the top left hand corner of 
Rp^q is positioned over the (^, j)-entry, M[i,j], of M, is given by summing all 
the entries on view: 

p-l q-l 

= ^^M[i + r,j+c], l<i<m-p+l, l<j<n-q+l. 

r=0 c=0 

Thus, we obtain an {m — p + 1) x {n — q + 1) matrix Rp^q(M) called the 
{p, q) -rectangular scan of M; when p and q are understood, we write R{M) — 
Rp^q{M), and speak more simply of the rectangular scan. (This terminology is 
a slight departure from that found in [6].) In the special case when R{M) has 
all entries equal, say k, we say that the matrix M is homogeneous of degree k, 
simply /c-homogeneous. 

Furthermore, we define an {m — p) x [n — q) matrix Xp^q[M) by setting, for 
1 < i < m — p, 1 < j < n — q: 

Xp,q{M) [i, 3] = M[i, j\ + M[i + p, j + g] - M[i + p, j] ~M[i,j+q]. 
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As usual, when p and q can be understood without ambiguity, we suppress 
them as subscripts. In the event that the matrix xi^) is a zero matrix, the 
matrix M is said to be smooth. Notice that the homogeneous matrices are 
properly included in the smooth matrices, as shown by the matrix M of Fig. 2, 
which is smooth, and non homogeneous. 
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Fig. 1. The matrix M and its corresponding matrices R2,3{M) and X2,3{M). 



Simplifying the rules of the game, thought the paper we will consider the 
matrix M (representing a physical structure) as a binary one; under this as- 
sumption, the rectangular scan R{M) turns out to be a positive matrix whose 
values are in the set {0, . . . ,p ■ q}, and the matrix x(M) turns out to have 
values in the set {—2, . . . , 2} (see Fig. 1). 

We conclude this introductory section with three observations which are direct 
consequences of the given definitions and to which we shall have frequent re- 
course in what follows. Since their proofs are a matter of simple computations, 
they are omitted. 

Lemma 1 // Mi and M2 are two m x n binary matrices, then 

R{Mi + M2) = R{Mi) + R{M2) and x{Mi + M2) = x{Mi) + x{M2) . 
Lemma 2 If M is a binary matrix, then 

XiARpAM)) = XpAM). 



Thus the rectangular scan R{M) of a binary matrix M already contains suf- 
ficient information to compute x(M) and so to decide whether M is smooth. 
Notice that, with a certain terminological inexactitude, we can also say, in 
the case where M is smooth, that R{M) is smooth (more precisely, R{M) is 
(1, l)-smooth, while M itself is (p, g)-smooth, as our more careful statement 
of the Lemma 2 makes clear). 

An appeal to symmetry and induction yields the following generalization of 
[6, Lemma 2.2]. 

Lemma 3 If M is a smooth matrix then, for any integers a and (3 such that 
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1 < i + ap < m and 1 < j + Pq < n, 

M[i, j] + M[i + ap, j + pq] = M[i + ap, j] + M[i, j + Pq] . 

Finally, we say that an entry M[i,j] of the matrix M is {p, q)-invariant if, for 
any integer a such that 1 < i + ap < m and 1 < j + aq < n, 

M[i + ap, j + aq] = M[i, j] . 

If all the entries of M are (p, g)-invariant, then M is said to be (p, g)-invariant. 



2 A Decomposition Theorem for Binary Smooth Matrices 

In this section we extend the studies about homogeneous matrices started in 
[6] to the class of smooth matrices: first we furnish a series of simple results 
which link smoothness and invariance, then we proceed along a path leading 
through a decomposition theorem for smooth matrices to their reconstruction. 

Lemma 4 If M is a smooth matrix, then each of its elements is {p, 0)-invariant 
or (0, q)-invariant. 

PROOF. Since M is smooth, for each 1 < i < m — p and 1 < j < n — g, it 
holds 

M[i, j] + M[i + p,j + q] = M[i + p, j] + M[i, j+q]. 
Let us consider the following three possibilities for the element M[i,j]: 

i) M[i,j] 7^ M[i + p,j]: by Lemma 3, for o; = 1 and for all /3 e Z such that 
I < j + Pq < n, it holds M[i,j + pq] = M[i,j] and M[i +p,j + Pq] = 
M[i + p,j], so M[i,j] is (0, g)-invariant. 

ii) M[i,j] 7^ M[i,j + q]: by reasoning similarly to i), we obtain that M[i,j] 
is {p, 0) -invariant. 

Hi) M[i,j] = M[i,j + q] = M[i +p,j]: if there exists cto G Z such that 
M[i + aQp,i] 7^ M[i,j], again reasoning as in i), we obtain that M[i,j] is 
(0, g)-invariant. 

On the other hand, if for all 1 < i + ap < m it holds that M[i + ap, j] — 
M[i,j], then M[i,j] is (p, 0)-invariant. 

Finally, ifm — p-|-l<i<m and n — q-\-l<j<n,&, similar reasoning leads 
again to the thesis. □ 

The reader can check that each entry of the smooth matrix M in Fig. 2 is 
(2, 0)-invariant (the highlighted ones) or (0, 3)-invariant. A first decomposition 
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result follows: 

Theorem 5 A matrix M is smooth if and only if it can be obtained by sum- 
ming up a {p,0) -invariant matrix Mi and a {0,q) -invariant matrix M2 such 
that they do not have two entries 1 in the same position. 



PROOF. (^) Let Ml and M2 contain the {p, 0)-invariant and the (0, q)- 
invariant elements of M, respectively. By Lemma 4, the thesis is achieved. 

(<^=) Since Mi is (p, 0)-invariant, then for each 1 < i < m — p, 1 < j < n — q 
it holds 

x{Mi)[i,j] = Mi[i,j] + Mi[i + p,j +q]- Mi[i + p,j] - Mi[i,j + q] ^ 
= Mi[i,j] + Mi[i,j + q]- Mi[i,j] - Mi[i,j + g] = 

so, by definition, Mi is smooth. The same result holds for M2 and, by Lemma 1, 
for M = Ml + M2. □ 

We can go further on by reformulating this last theorem in terms of the rect- 
angular scans of the matrices Mi and M2 : 

Lemma 6 The following statements hold: 

i) if M is if), q) -invariant, then R{M) has constant rows; 

ii) if M is {p,0) -invariant, then R{M) has constant columns. 

PROOF, i) For each l<i<m— p-\-l and I < j < n — q, we prove that 
R{M)[i,j]^R{M)[i,j + l]: 

R{M) [z, J + 1] = E^;=o E^=o M[t + r,j + l + c] = 

= E?=o E?=! M[i + r,j + c]+ E?=o M[i + r, i + g] = 
since M is (0, g) -invariant 

= E E M[i + r,j + c\ + M[z + r, j] = R{M) [i, j] . 
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ii) The proof is similar to i). □ 



After observing that each matrix having constant rows or columns is smooth, 
a direct consequence of Theorem 5 and Lemma 6 is the following: 

Theorem 7 A binary matrix M is smooth if and only if R{M) can be de- 
composed into two matrices Ry. and Rc having constant rows and columns, 
respectively. 

Fig. 2 shows that the converse of the two statements of Lemma 6 does not 
hold in general. However, we can prove the following weaker version: 

i\ r 1 

}': 1 1 H? 3': 3 J 



\ \ \1 ' 1 l: 1 1 

1\ 1\ \ 1 2 2\ 2 2 

Fig. 2. A non invariant matrix M whose {2, 3) -rectangular scan has constant rows. 
Proposition 8 Let M be a binary matrix. The following statements hold: 

i) if R{M) has constant columns, then there exists a {p,0) -invariant matrix 
Ml such that R{M) = R{Mi); 

ii) if R{M) has constant rows, then there exists a {0,q) -invariant matrix M2 
such that R{M) = R{M2). 

PROOF, i) We define the matrix Mi as follows: the first p rows of Mi are 
equal to those of M, and the other entries of Mi are set according to the 
desired (p, 0)-invariance. It is easy to verify that R{Mi) = R{M). 

ii) A definition of M2 similar to that in i) for Mi can be easily given. □ 

2. 1 Solving Reconstruction (A, p, q) for Smooth Matrices 

A first approach to the general reconstruction problem consists in the defini- 
tion of the following algorithm which suites only for a binary smooth matrix 
whose (p, g)-rectangular scan has constant rows: 

RecConstRows(A, p, q) 

Input: an integer matrix A of dimension m' x n', having constant rows, and 
two integers p and q. 
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Output: a (0, g)-invariant matrix M, of dimension m x n, where m = ni' + 
p — 1 and n = n' + q — 1, having A as (p, q')-rectangular scan, if it exists, 
else return FAILURE. 

Procedure: 

Step 1: create the mxn matrix M and the vector PEnt (storing the Partial 
number of Entries 1 in each row of M) of dimension m, to support the 
computation. Initialize the entries both of M and PEnt to 0. 
For each row 1 < i < m' — 1, 
Step 1.1: if A[i, 1] < A[i + 1, 1] then 

• M[i + p,l] = --- = M[i+ p, A[i + 1,1]- A[i, 1] + PEnt\i]] = 1; 

• PEnt[i + p] = A[i + 1, 1] - A[i, 1] + PEnt[i]. 
If PEnt[i +p]> q then FAILURE. 

Step 1.2: if A[i, 1] > A[i + 1, 1] and PEnt\i] > A[i, 1] - A[i + 1, 1] then 

• M[i+p,l] = ---^ M[i+ p, A[i + 1,1]- A[i, 1] + PEnt\i]] = 1; 

• PEnt[i + p] = A[i + 1,1]- A[i, 1] + PEnt[i]. 

Step 1.3: if A[i, 1] > A[i + 1, 1] and PEnt\i] < A[i, 1] - A[i + 1, 1] then 

• k = A[i, 1] - A[i + 1,1]- PEnt\i]] 

• for each i' < i, i' = {i)modp 

. M[i', PEnt[i'] + 1] = ■ ■ ■ = M[i', PEnt[i'] + k] = 1; 

. PEnt[i'] = PEnt[i'] + k; 

. if PEnt[i'] > q then FAILURE. 

Step 2: let k = A[l, 1] - PEnt[l] PEnt\p\. 

For each 1 <k' <k, search one of the upper leftmost pxq positions of M, 
say such that, for each i' = {i)modp, \ <i' < m, it holds M[i',j] — 0. 

If such a position does not exist then FAILURE, 
else set all the entries M[i',j] to the value 1, and increase k' by one. 

Step 3: complete the entries of M according to the (0, g)-invariance con- 
straint, and return M as OUTPUT. 

As regard the correctness of this reconstruction algorithm, it rehes on the 
analysis of what stored in M after Step 1: at that stage, in fact, the entries in 
first column of its rectangular scan R[M) differ from those of A by the same 
constant value, without overcoming. 

The formal counterpart of what sketched above is in the following lemmas: 

Lemma 9 After performing Step 1 o/RecConstRows(A,p, q), for each 1 < 
i < m' , it holds 

R{M)[i, 1] - R{M)[i + 1,1] = A[i, 1] - A[i + 1,1]. 
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PROOF. Let us first inspect the entries placed in the rows i and i + p of M 
during Step 1 of RecConstRows(A,p, q), for a generic index 1 < z < m': 

if A[i, 1] < A[i + 1,1], then in row z + p of M are added A[i + 1,1]- A[i, 1] + 
PEnt[i] entries 1. Since row i of M contains PEnt[i] entries 1, then, at that 
step, it holds 

R{M)[t,l] - R{M)[i + 1] = A[i,l] - A[z + (1) 

if A[i, 1] > A[i + 1,1], and PEnt[i] > A[i, 1] - A[i + 1, 1], then in row i + p 
of M are added PEnt[i] — A[i, 1] + A[i + 1,1] entries 1, so equation (1) still 

holds; 

if A[i, 1] > A[i + 1, 1], and PEnt[i\ < A[i, 1] - A[i + 1,1], then in row i of M 
arc added A[i, 1] — A[i + 1,1]— PEnt[i] entries 1 in addition to the PEnt[i] 
ones already present, so equation (1) is again satisfied. 

For each i < i' < i + p, Step 1 eventually changes some entries from row i + 1 
to row i + p — 1 of M. So, both R{M)[i, 1] and R{M)[i + 1, 1] increase their 
value of the same amount, without compromising the validity of equation (1). 

Finally, if i' > i + p, then Step 1.3 may modify the vahics of R{M)[i, 1] and 
R{M)[i + 1, 1], but again of the same amount, since the (eventually) added 
entries 1 respect the {p, 0)-invariance in the rows of index less than i', so 
equation (1) definitively holds, and we obtain the thesis. □ 

Lemma 10 Let us consider the vector PEnt as updated at the end of Step 1 
o/ RecConstRows(A,p, g). For each matrix M such that R{M) = A, and 
for each < i < m, it holds 

M[i,l] H hM[i,p] > PEnt\i]. 



PROOF. By contradiction, we assume that there exists an index 1 < i < m 
and a matrix M' such that 

M[i, 1] H h M[i,p\ + PEnt\i], 

with R[M) = A and A; > 0. By Lemma 9, the same equation holds for each 
row i' — {i)modp. 

Let Zo be the first index such that 

• iQ — {i)modp; 

• for each i' — {i)modp, it holds PEnt[io\ < PEnt[i']. 
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If < m — p, then the minimality of the value of PEnt[io\ assures that 
PEntlio] = before Step 1 reached row io, and, consequently, that A[iQ, 1] < 
4^0 + 1,1]. 

As soon as Step 1.1 reaches the row index iq, it eventually increases the value 
PEnt[iQ + p], leaving unchanged that of PEnt[iQ]. Again the minimality of 
PEnt[io] assures that no changes will be performed to the value of PEnt[io] 
till the end of Step 1. 

Hence, the assumption M[io, 1] + • • • + M[io,p] + k — PEnt[io\, with A; > 0, 
generates a contradiction. 

li io > m — p, then a similar argument holds, and so we get the thesis. □ 

Now, also Step 2 of RecConstRows {A,p, q) can be better understood: the 
k elements which change their value to 1 and which are added in the upper 
leftmost p X q positions of M, fill the gap among R{M)[i. 1] and 74[i,l], so 
that the output matrix M has the desired property R{M) = A. 

Corollary 11 Each row i of a matrix M having A as rectangular scan con- 
tains at least PEnt[i] elements which are {p,0) -invariant, and not (0, g)- 
invariant. 

A procedure which reconstructs a smooth matrix whose {p, g)-rectangular scan 
A has constant columns, say RecConstCols(A,p, g), can be easily inferred 
from RecConstRows, so, in the sequel, we will consider it as already defined. 

Prom Lemmas 9 and 10, it is straightforward that 

Theorem 12 The problem Reconstruction{A, p, q) can be solved in 0{mn), 
when A has constant rows or constant columns. 

Example 13 Let us follow the computation RecConstRows (^,3,4), with A de- 
picted in Fig. 3. 
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Fig. 3. The matrix A of Example 13. 

Step 1 .• the matrix M is created, and its first four columns, together with the vector 
PEnt, are modified as shown in Fig 4- More precisely, 

A[l, l]+2 = ^[1, 2] requires Step 1.1 to place two entries 1 in (the leftmost positions 
of) row 4, Fig 4, (a); 
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A[l,2] = A[l,3] + 1 requires Step 1.3 to place one entry 1 in row 2, Fig 4, (b); 

A[l,3] + 2 = ^[1,4] requires Step 1.1 to place two entries 1 in row 6, Fig 4, (c); 

A[l, 1] = A[l, 2] + 3 requires Step 1.3 to add one entry 1 both in row 1 and in row 
A, Fig 4, (d). 



I I 1:1::: 11::': 1:11 



I I I I 



PEnt: (0,0,0,2,0,0,0) PEnt: (0,1,0,2,0,0,0) PEnt: (0,1,0,2,0,2,0) PEnt: (1,1,0,3,0,2,0) 
(a) (b) (c) (d) 

Fig. 4. Step 1 o/RecConstRows(A,3,4). 

Step 2 places the remaining k = A[l, 1] — PEnt[l] — PEnt[2] — PEnt[3] = 3 entries 
1 in the upper leftmost px q suhmatrix of M , and propagates them according to the 
{3,0)-invariance, paying attention that no collisions occur (see Fig. 5 (a)). 

Step 3 completes M according to the {Q,A)-invariance, giving the final solution de- 
picted in Fig. 5 (6). 
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(a) (b) 
Fig. 5. Step 2 and Step 3 o/ RecConstRows(^, 3, 4). 

The general reconstruction algorithm 

Theorem 7 allows one to foresee the use of the procedures RecConstRows 
and RecConstCols to solve Reconstruction g), when A is (1, l)-smooth: 
the algorithm at first will split the matrix A into two parts having constant 
rows and columns, respectively, then it will apply to each of them the ap- 
propriate reconstruction procedure, and finally it will merge the two outputs. 
Performing the merging stage a conflict occurs when the same position in two 
output matrices has value 1. To prevent it small refinements to the outputs of 
RecConstRows and RecConstCols will be required. 
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So, let us start by showing in the next lemma a qiiick way of finding all 
the possible decompositions of a (1, l)-smooth matrix into two parts having 
constant rows and columns, respectively 

Lemma 14 Let A be a m x n integer matrix. If A is (1, l)-smooth, then it 

admits k + 1 different decompositions into two matrices having constant rows 
and columns, with k being the minimum among all the elements of A. 

PROOF. The thesis is achieved by defining a procedure which gives as output 
a complete list of couples of matrices (A*,^*), with < t < A;, each of them 
representing a decomposition of A into two parts having constant rows and 
columns, and successively, by proving its correctness: 

Decompose (A) 
Input: an integer m x n matrix A. 

Output: a sequence of different couples of matrices {A^, A^), . . . , {A^, A^), 
with k being the minimum element of A, such that, for each 1 < t < k, A^ 
has constant rows, A* has constant columns, and Al + Al — A. If such a 
sequence does not exist, then return FAILURE. 

Procedure: 

Step 1: initialize all the elements of two m x n matrices Ac and Ar to the 
value 0. Let k be the minimum among the entries of A. From each element 
of A, subtract the value k and store the result in A^, 

Step 2: for each 1 < i < m 
Step 2.1: compute 

ki = minj{Ac[i,j] : 1 < j < n}; 

Step 2.2: subtract the value ki from each element of A^, 
Step 2.3: set all the elements of row i of A^ to the value kf, 
Step 3: if the matrix Ac has not constant columns then FAILURE 
else for each <t < k, create the matrices A* and ^4* such that 

j] = Ar[t, j] + t and A% j] = Ac[i, j] + k- t, 

with 1 < i < m and I < j < n. 

Give the sequence (A^, ^°), . . . , (A^, Aj) as OUTPUT. 

Example 15 shows a run of the algorithm. By construction, each couple (A* , A*) 
is a decomposition of A, and furthermore, Al has constant rows. 

What remains to prove is that the matrix Ac updated at the end of Step 2 has 
constant columns (and, consequently, the same hold for all the matrices A*). 
Let us denote by rj the common value of the elements of the z-th row of A^., and 
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let us proceed by contradiction, assuming that Ac has not constant columns. 
Since A is the sum of a column constant and a row constant matrix, and for all 
1 <i < m — p+l and 1 < j < n — q + 1, it holds j] = ri + Ac[i,j] + k, then 
Ac is also the sum of a column constant matrix, and a row constant matrix, 
this last having at least one row, say io, whose elements have value k'^^ ^ 0. 

This situation generates an absurd, since ki^^ computed in Step 2.1 turns out 
no longer to be the minimum of row io in Ac, updated to that step. 

Since a matrix having constant rows (resp. columns) cannot be obtained as 
sum of a matrix having constant rows and a matrix having constant columns 
unless the latter is a constant matrix, then the k + 1 decompositions listed by 
the algorithm are all the possible ones. □ 

Example 15 Let us follow the steps of the procedure Decompose(A), with the 
matrix A depicted in Fig. 6. 
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Fig. 6. The {l,l)-smooth matrix A of Example 15. 

Step 1: we subtract from all the elements of A, the value k = 2, i.e. its minimum 
element, and we store the obtained result in the matrix Ac. 

Step 2: for each 1 < i < m—p+l, we find the minim.um, value hi am,ong the 
elements of row i of Ac (Step 2.1), we subtract it from all these elements (Step 2.2), 
and finally, we set the elements in row i of A^ to the value ki (Step 2.2>). In our 
case, the minimums are k\ =2, k2 = 1, k^ = Q and A;4 = 1. 
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Fig. 7. The three decompositions of the matrix A. 

Step 3: the matrix Ac updated at the end of Step 2 has constant columns, so the 
three different decompositions of A can be computed and listed. 



The output is depicted in Fig. 7. 
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Now wc arc finally able to define a general procedure which solves the problem 
Reconstruction {A,p,q), when A is (1, l)-smooth: 



RecSmooth {A,p,q) 

Input: an integer (1, l)-smooth matrix A of dimension m' x n' and two 
integers p and q. 

Output: a binary matrix M of dimension m x n, with m = m' + p — 1 and 
n = n' + q — 1, having A as (p, g)-rectangular scan, if it exists, else return 
FAILURE. 

Procedure: 

Step 1: run Decompose(^), and let ^4°), . . . , (A^, A^) be its output. 
Set t = 0; 

Step 1.1: run Step 1 of RecConstRows g). 

Let PEntrow = PEnt and k,.ow = k, with PEnt and updated at the 
end of the step. Define PEnt=p to be the vector having p elements, and 
such that: 

PEnt=p[i] — 'max{PEntrow['i''] '■ i' — {i)modp}, 

with 1 < i < p, and 1 < i' < m. 
Step 1.2: run Step 1 of RecConstCols(A*,p, g). 

Let PEntcoi = PEnt and kcoi = k, with PEnt and k updated at the 
end of the step, and let PEnt=q be the vector having q entries, and such 
that: 

PEnt=q[j] = max{PEntcoi[j'] : f = {j)modg}, 
with 1 < J < g, and 1 < f < n; 

Step 1.3: among all the possible 2?xg matrices whose entries are in {P, Q, 1, 0}, 

choose one, say W, such that: 

i) the number of entries Q in its i-th row is PEnt=p[i]. 

ii) the number of entries P in its j-th column is PEnt=q[j]; 
Hi) the number of entries 1 is krow + kcoi- 

liW does not exist and t ^ k then set t = t + 1, and return to Step 1.1. 
If W docs not exist and t = A; then FAILURE; 

Step 2: create the m x n matrix M, and initialize its entries as follows: 

Step 2.1: for each < i < m and for each < j < g, 

if PEntrow[i] 7^ and = Q, with i' = {i)modp then set both 

= Q and PEntrow[i] — PEntrow[i] — 1; 
Step 2.2: for each < j < n and for each < i < p, 

if PEntcoiij] 7^ — P, with / = {j)modq then set both 

M[i,j] = P and PEntcoiij] = PEntcoi[j] - 1; 
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Step 2.3: for each < z < p and for each < j < g, 
if W[iJ] = 1, then set M[iJ] = 1; 

Step 2.4: fill the matrix M imposing the (0, g)-invariance of its entries Q 
and 1, and the (p, 0)-invariance of its entries P and 1; 

Step 3: change the values P and Q to 1, and set the remaining entries of M 
with the value 0; finally, give M as output. 

Theorem 16 The problem Reconstruction {A,p, q), with A being (1, \)-smooth, 
admits a solution if and only z/RecSmooth [A,p,q) does not return FAIL- 
URE. 



PROOF. (^) Let M be a solution of Reconstruction g), and let us 

assume that M = Mi + M2, with Mi and M2 being (0, g)-invariant and {p, 0)- 
invariant, respectively. 

Let (A*, A*) be one of the decomposition of A such that R{Mi) — A* and 
R{M,) = Al 

Lemma 10 implies that, for each 1 < i < m, the value PEntrowii] indicates the 
minimum number of elements of M which are (0, g)-invariant and not {p, 0)- 
invariant, and which lie in the first q columns of the solution; a symmetrical 
property holds for PEntcoi- 

Let us construct a. p x q matrix W' as follows: 

- for each 1 < i < m and 1 < j < g, if -M[z, j'] = 1 is (0, g)-invariant and not 
(p, 0)-invariant, then set j] = Q, with i' — {i)modp\ 

- for each 1 < i <p and 1 < i < n, if M[i,j] = 1 is (p, 0)-invariant and not 
(0, g) -invariant, then set W'[i',j] = P, with j' = {j)rnodg; 

- for each 1 < i < p and 1 < j < g, if M[i,j] = 1 is (p, 0)-invariant and 
(0, g) -invariant, then set j] = 1. 

Obviously, by definition of invariance, in W there are no positions which are 

first set to a value and then modified to another. So, the existence of matrix 
W implies that of a matrix W satisfying the constraints imposed in Step 1.3. 

("4=) Immediate. □ 

Example 17 Let us describe a run of RecSmooth {A, 3, 3) starting from 
the decomposition of A into the couple of matrices depicted in Fig. 8. 
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Fig. 8. The decomposition of the rectangular scan A used in Example 17. 
Step 1.1 produces the vectors 

P£;nt^<^ = (0,0,2,0,1,0,0) and PEnt=p = (0,1,2), 

while Step 1.2 produces the vectors 

P£;ntcoZ = (0,1,1,0,0,0,0,0,1,1) and PEnt=q ^ (1,1,1). 
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Fig. 9. The matrices created in successive steps o/RecSmooth (^,3,3). 

Among the matrices which are compatible with the requirements of Step 1.3 we 
choose that depicted in Fig. 9, (a). 

Steps 2.1, 2.2 and 2.3 produce the matrix M in Fig. 9, (h) (notice that the 

entries P are (p,0) -invariant, while the entries Q are (0,q)- invariant as de- 
sired). Finally, Step 2.4 produces the matrix in Fig. 9, (c), and, consequently. 
Step 3 the output. 

The following theorem holds: 

Theorem 18 The computational complexity o/REcSMOOTH(yl.,p, g) is poly- 
nomial in m and n. 



PROOF. We obtain the thesis by analyzing the complexity of each step of 
RecSmooth: 

Step 1: the procedure DECOMPOSE clearly acts in 0(mn) time (remind that 

k < p ■ q is the minimum among the elements of ^4) ; 
Step 1.1 and Step 1.2 are repeated at most k times, and, each time, they 

ask for a run of RecConstRows and of RecConstCols which are both 
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performed in 0{mn). The compTitation of PEnt=p and PEnt=q does not 

increase the complexity of these two steps. 
Step 1.3 is carried on in constant time with respect to m and n. 
Steps 2 and 3 require 0{m n) to create matrix M. 

Hence, the total amount of time is 0(m n). □ 

Remark 19 We are aware that Step 1.3 o/ RecSmooth {A,p,q), i.e. the 
search of the matrix W , can he carried on in a smarter way, hut this will hring 

no effective contribution to the decreasing of the computational complexity of 
the reconstruction, and, on the other hand, it will add new lemmas and proofs 
to the current section. 

The first part of the paper devoted to the analysis and the reconstruction of 
smooth matrices is now completed. 



3 Solving Reconstruction (v4,p, g): final challenge 

This last section concerns the matrices which are not smooth: in particular, 
for each non smooth matrix M, we consider the matrix x(M) and we define 
a polynomial time algorithm which lists all the matrices consistent with it. 
Finally we will integrate it with the algorithm for reconstructing a smooth 
matrix defined in the previous section, and we will achieve the solution of the 
general reconstruction problem. 

Unfortunately, the definitions introduced up to now are not specific enough 
to describe these further studies, and a final effort is required to the reader: 
what follows has the appearance of a stand alone part inside this section, but 
the feeling of a final possible usage will never be frustrated. 

Hence, let a and b be two indexes such that l<a<p,l<b<q, and A be 
an integer mxn matrix. We define the (a, b)-subgrid of A to be the submatrix 

SiA)a,b[hj] = A[a + ii-l)p,b+ (j - 1) q] 
with 1 < a + (i — l)p < m and 1 < b + {j — l)q < n (see Fig. 10). 
If we consider again a binary matrix M, by definition it holds that 

x{M)[a + (i - 1) 6 + (j - 1) g] - S{x{M)\,[i,j] = 

= S{M),,,[i,j] + S{M),,,[i + 1, j + 1] - S{M\,,[i + 1, j] - S{M),,,[i,j + 1]. 
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The binary matrix V of dimension m x n is said to be a valuation of S{x{M))a,b 
if, for each l<i<m, l<j<n, 

- if i ^ {a)modp and j ^ {h)modq then j] = 0; 

- S{x{M))a,b = S{xiy))a,b (see Fig. 10). 

The notion of valuation extends to the whole matrix x(M) as the union of the 
valuations of all its subgrids. 

Proposition 20 Let S{x{.^))a,h and S{x{M))a',b' be two subgrids whose val- 
uations are V and V , respectively. If a ^ a' orb^ b' , then for each 1 < i < m 
and 1 < j < n, V[i,j] = 1 implies V'[i,j] = . 

7 ; ; : 7 ; : 1 
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5(M) 



Si X (M)) 



SiV) 



Fig. 10. T/ie subgrids of the matrices M and xi^) with respect to the position 
(2,2). Matrix V is one of the possible valuations of S{x{M))2^2- 

Lemma 21 Let V be a valuation of S{x{^))a,h, and let Zq be a row [column] 
of S{V)a,b having all the elements equal to 1. The matrix V such that S{y')a,b 
is equal to S{y)afi except in the elements of the row [column] iq which are all 
set to 0, is again a valuation of S{x{.^))a,b- 

If V and V' are two valuations as in Lemma 21, then we say that the valuation 
V is greater than the valuation V' . This relation can be easily extended to a 
finite partial order on the valuations of the subgrids of x{^)- 



Lemma 22 Let 1 < i < m — p and I < j < n — q. //x(M)[i, j] 
M[i,j] = M[i + p,j + q] = l, and M[i + pj] = M[i,j + q]=0. 



2, then 



The proof is immediate. A symmetric result holds if x{M)[i,j] has value —2. 



The following lemma turns out to be crucial in this section. A deeper analysis 
of what it states could furnish material for further studies: 



Lemma 23 Given a binary matrix M, for each couple of integers 1 < a < 
1 < 6 < g, the number of minimal elements in the partial ordering of the 
valuations of S{x{.^))a,b is polynomial with respect to the dimensions m and 
n of M. Furthermore, each minimal element can be reconstructed in polynomial 
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Fig. 11. The four valuations of the point pi = {2, 2) of value 1. 
time with respect to them. 



PROOF. 



Let S{x{M))a,b have dimension m' x n' . We order the (positions of the) non 
zero elements of S{x{^))a,b according to the numbering of its columns (from 
left to right), and, in the same column, according to the numbering of its rows 
(from up to bottom), and let pi, . . . , pt be the obtained sequence. We prove the 
thesis by induction on the number t of elements of the sequence, i.e. we prove 
that the addition of new nonzero elements in S{x{^))a,b does not increase 
"too much" the number of its possible minimal valuations. 

We first observe that, by Remark 22, the presence of entries 2 or —2 in 
S{x{^))a,b does not increase the number of minimal valuations, so we are 
allowed to focus our attention exactly on the elements of value 1 or —1. As 
one can expect, the symmetry of the two cases allows us to show the details 
of one (in particular when the element to add has value 1), and let the reader 
infer the other: 

Base t = 1: if = (z, j) and S{x{.^))a,b[h j] = 1, then the four possible 
valuations of S{x{M))a,b are depicted in Fig. 11. Among them, only S(yi)a,b 
and S{V2)a,b are minimal: they can be reached both from S(ys)a,b and from 
S{V4)a,b by deleting the rows or the columns entirely filled with entries 1, 
as stated in Lemma 21. 

Step t — > t + 1: let S{x{M))a,b have a sequence pi, . . . ,pt+i of nonzero 
points, with pt+i = and S (xi^)) a,b[i , j] = 1- Let V he a valuation 

of the first pi, . . .pt points in S{x{M))a,b- It is straightforward that, for all 
1 < < m' + l and J < / < n' + l it holds S{V)a,b[i' , J + ^] = S{V)a,b[i' , f]- 
Hereafter, we show all the possible ways of extending V to the valuation 
V which includes the point Pt+i- Some pictures are supplied in order to 
make the different cases transparent. 
Let us call 

0-row: a row of S(y)a,b whose elements have all value 0; 

*l-row: a row of S(y)a,b whose element in column j + I has value 1; 
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*0-row: a row of S{V)a,b which is neither 0-row nor 1-row. 
We examine all the possible configurations of S{V)a,b, and for each of 
them we indicate the desired extension to the minimal valuation S{V')a,b- 

i) all the rows from 1 to i are 0-rows or *l-rows (see Fig. 12, (a)). 

We define the valuation S{V')a,b as follows: for each 1 < i' < i, 
if row i' is a 0-row, then change from to 1 the value of each entry of 
S{V)a,b in position with 1 < j' < j, so that it becomes a *0-row; 

if row i' is a * 1-row, then change from 1 to the value of each entry of 
S{V)a,b in position {i',j'), with j + 1 < j' < n'+q—l, so that it becomes 
a *0-row or a 0-row. If a 0-row is created, then discard the obtained 
valuation S{V')a,b, since it has been already obtained in a previous step 
(easy check); 

ii) there exists a *0-row ^', with 1 < i' < i (see Fig. 12, (6)). 

No changes in the first i rows of S{V)a,b allow the insertion of the new 
point pt+i, 

in) all the rows from i-|-ltom'-|-p — 1 are 0-rows or *0-rows (see Fig. 12, 
(c)). 

We define the valuation S{V')a,b as follows: for each i + 1 < i' < 
m' + p — 1, change from to 1 the value of each entry of S{V)a,b in 
position with j + 1 < j' < n' + q — 1, so that it becomes a *1- 

row. Discard such a valuation if a row having all the entries equal to 
1 has eventually been created, in order to maintain minimality (remind 
Lemma 21); 

iv) there exists a *l-row i', with i + l<i'<n' + q — 1 (see Fig. 12, (d)). 
No changes in the last (m' +p—l)—i rows of V allow the insertion of 

the new point pt+i- 

Remark: If the point pt+i has value —1, then a further check is needed in 
the analogous of case i): it may happen that there exists a *0-row i', with 
1 < i' < i, which is turned into a *l-row after the addition of Pt+i, so that 
a non minimal configuration is created. 

The four configurations above described arc exhaustive with respect to 
the addition of the single point pt+i of value 1 to the valuation S{V)a,b- 
However, has not yet been considered: it appears when two points 

Pt^i and pt+2 are added to S{V)a,b, under the assumption that they have 
different value, and they lie in the same column. 

Step t t+2: let us assume that pt+i = {i,j),pt+2 = S{x{M)a,b[i, j] = 

1, and S{x{M))a,b[i',j] = -'i- 
If it holds that: 

v) all the rows from i to i' are not *l-rows of S{V)a,b, and there exists 
a *l-row with index greater than i' (this last condition prevent S{V)a,b 
from being extended to S{V')a,b by means of Hi)). 

We define the valuation V as follows: 

for each i < io ^ i', change from to 1 the value of each entry of 
S{V)a,b in position (io, /), with j -|- 1 < / < n' -|- g — 1, so that it becomes 
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Fig. 12. Examples of the possible ways of extending a valuation V when adding 
point pt+i- In cases (b) and (d), the boldface rows prevent V from being extended. 



a *l-row. 



In the sequel, when we mention the above described cases i) — v), we intent 
to include also their symmetrical counterparts. It is immediate to check that 

- cases i) — v) extend S{V)a,b by adding the desired point (or points); 

- each extension of S{V)a,b is minimal, since no rows or columns completely 
filled with entries 1 are added; 

- all the minimal valuations for the sequences pi, ■ ■ ■ ,Pt+i or pi, . . . ,pt+i,pt+2 
are obtained by means of i) — v). 

So, what remains to prove is that the number of the different minimal valu- 
ations for a given matrix S{x{M))a,b is polynomial in its dimensions m' and 
n' (and consequently in the dimensions m and n of M). We achieve this aim 
by showing that the number of valuations S{V)a,b which admit more than a 
single minimal extension is bounded by m'. Some properties are needed: 

Proposition 24 Each valuation S{V)a,b admits at most two different mini- 
mal extensions both when adding a single point pt+i (see Fig. 13), and when 
adding two points pt+i o,ndpt+2, under the assumptions ofv). 

Proposition 25 // the valuation S{V)a,b does not contain any 0-row, then it 
admits at most one minimal extension via i) — v) (see Fig. 13, valuation U). 

Proposition 26 The 0-rows of each valuation which extends S{V)a,b are a 
subset of those of S{V)a,b- Furthermore, if S(y)a,b extends in two minimal 
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Fig. 13. The valuations S{V)a,b Oind S{U)a,b extend in two different valuations. 
Since S{U)a,b does not contain any 0-row, one of its extensions (i.e. S{U2)a,b) is 
not minimal. 

ways, then the two extensions do not share any 0-row. 

Proposition 27 Two minimal valuations of S{x{M))a,b do not share any 0- 
row. 



The proofs of these properties directly follow from the definitions of i) — v). 

Hence, Property 27 assures that each matrix S{x{M))a,b has at most m' dif- 
ferent minimal valuations containing 0-rows. From Property 25, if we add one 
or two new points to them, then at most m new minimal valuations may arise. 
As a neat consequence, we obtain that the number of minimal valuations of a 
given matrix S'(x(M))a f, is polynomial in m' and n', and so it is the complexity 
of their reconstruction. □ 

Example 28 Let us find all the minimal valuations of the matrix S{x{M))a,b 
depicted in Fig. I4 



1 
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-1 



Fig. 14. The matrix S{x{M))a,b of Example 28. 

We proceed from the leftmost entry of S{x{M))a^b different from 0, till the rightmost 
one, and we construct, step by step, all the possible minimal valuations, as described 
in i) — v). 

The computation is represented in Fig. 15, by using a tree whose root is the matrix 
having all the entries equal to 0, and which represents the (a, b) subgrid of the 
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valuation of a 0-homogeneous matrix. The nodes at level k are all the possible (a, b) 
subgrids of the valuations of the first k entries different from of S{x{^))a,b- 



On each matrix, the highlighted cells refer to the correspondent entry 1 or —1 of 
5(x(Af))a,b which is being considered. 
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Fig. 15. The computation of the minimal valuations of S{x{M))a^h- 

It is easy to check that any further addition of 1 or —1 entries in S{xiM))a^b does 
not increase the number of the minimal valuations. 

The following variant of the procedure RecSmooth will be used in the final 
reconstruction algorithm; the details of the procedure which differ from the 
original ones are given: 



RecSmoothAll(v4,]?, q) 

Input: an integer matrix A and two integers p and q. 
Output: a (eventually void) sequence of m x n matrices Mi, 



(1) 



elements in {Ip 

Procedure: 

Step 1: ... 

Step 1.1: ... 



(fccoi) 1 (1) 



1 (^roii.') 



Mk having 



, P, Q, 0}, with 1 < n < p X q. 
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Step 1.2: ... 

Step 1.3: list all the possible p x q matrices such that, for each of them 

i) the number of the entries Q in its i-th row is PEnt=p[i]. 

ii) the number of the entries P in its j-th column is PEnt=q[j]; 

in) at least one occurrence of each entry in Ip-* , . . . , 1^=°'^ and at least 
one occurrence of each entry in 1^\ . . . , 1^''°"") is present. Furthermore, 
all the entries with pedex P [resp. Q], and having the same index must 
lie in the same column [resp. row]. 

Ii t ^ k then set t = t + 1, and return to Step 1.1. 

Step 2: let Wi, . . . Wk be the output list of Step 1. For each I < i < K, use 
matrix Wi to create the m x n matrix Mj whose entries are initialized as 
follows: 

Step 2.1: ... 

Step 2.2: ... 

Step 2.3: for each < i' < p and < j < q, 

iiWi[i',j] ^ {P,Q}, then set M^'J] = W,[i',j]]; 
Step 2.4: fill the matrix M imposing the (0, g)-invariance of its entries 
Q, I?,..., l''p^°'\ Iq\..., ljJ'-°"'\ and the (p, 0)-invariance of its entries 

P 1 (1) -\{kcol) 1 (1) -]{kro-w). 

I-, Ip , . . . , Ip , 1q , . . . , 1q , 

Step 3: return the sequence Mi, . . . Mk as output. 

This variant of RecSmooth inherits its 0{mn) computational complexity. 

As a final observation, one may wonder the meaning of the indexed entries 
Ip and 1q inside each matrix: the elements Ip [resp. 1q] having the same 
index mark the positions where a set of entries whose rectangular scan is 1- 
homogeneous, can be placed. The choice of all the possible sets of positions 
for the placement of the krow + A;corhomogeneous part of A constitutes a key 
point in the definition of the final reconstruction algorithm which follows: 

Reconstruction(A,P, q) 
Input: an integer matrix A and two integers p and q. 

Output: an m X n binary matrix M having A as {p,q) rectangular scan, if 
it exists, else return FAILURE. 

Procedure: 

Step 1: for each 1 < a < p and 1 < b < q, compute the sequence of minimal 
valuations of S{xi,i{A))a,b', 

Step 2: sum in all possible ways an element from each sequence of valuations 
computed in Step 1, and let Mi, . . . My be the obtained sequence of binary 
matrices; 
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Step 3: for each 1 < t < v, 

Step 3.1: compute the matrix At = A — R{Mt); 

Step 3.2: run REcSMOOTHALL(At,p, q), and let M{, . . . , be its output. 

Set t' = 1; 

Step 3.3: until t' < K, compute a matrix M by merging the matrix Mt 
and the matrix M/, as follows: initialize M — Mt; 
for each 1 < i < m, 1 < j < n 
if M;,[i,i] e {P,Q}, then 

if M[i,j] — 1, then set t' — t' -\- 1 and return to Step 3.3, else 
M[i,j] = l; 
if M;,[z, j] = l^p''\ with 1 < no < kcou then 

if M[i' ,j] = 0, for each position with i' = {i)modp, then set 

M[i', j] = 1, and change to the remaining entries in column of M^, 
having value lp°^ else set Mt,[i',j] = 0, and, if no other elements lp°^ 
are in column j, set t' = t' + 1 and return to Step 3.3; 
if Ml,[i,j] = 1q^°\ with 1 < mo < Kow, then 

if M[i,j'] = 0, for each position (i,/), with / = {j)modq, then set 
M[i,/] = 1, and change to the remaining entries in row i of M/, 
having value else set M/,[i,/] = 0, and, if no other elements 
are in row i, set t' = t' + 1 and return to Step 3.3; 
Return matrix M as output; 
Step 4: return FAILURE. 

The correctness of the procedure is straightforward, since we create all the 
possible minimal valuations for the entries of A which prevent it from being 
smooth, and successively, we merge them with all the possible solutions for its 
remaining smooth entries. 

Theorem 29 The problem Reconstruction {A,p, q) admits a solution, if and 
only if the algorithm Reconstruction {A, p, q) finds it. 

However, one can ask whether such a search always produces an output in an 
amount of time which is polynomial in the dimensions m and n of the solution. 
The answer is given in the proof of the following 

Theorem 30 The computational complexity o/Reconstruction(A,p, 5) is 
polynomial in the dimension m x n of the solution. 



PROOF. The complexity of the algorithm can be computed as the sum of 
the complexities of its steps, in particular: 

Step 1: Lemma 23 assures that the computation of all the valuations of 
S{xi,i{^))a,b can be performed in polynomial time with respect to m and 
n. 
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Step 2: the procedure asks for summing in all possible ways an element from 
each of the p ■ q sequences of valuations created in Step 1. Since each sum 
is performed in 0{m n), then the total complexity remains polynomially 
bounded by m and n. 

Step 3: the computation of the matrix At, the polynomial procedure Rec- 
SmoothAll and the merging process of Mt with Mj., are performed a 
polynomial number of times, without increasing the total complexity of the 
algorithm. □ 
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